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The "close hmit," a method based on perturbations of Schwarzschild spacetime, has proved to be 
a very useful tool for finding approximate solutions to models of black hole collisions. Calculations 
carried out with second order perturbation theory have been shown to give the limits of applicability 
of the method without the need for comparison with numerical relativity results. Those second order 
calculations have been carried out in a fixed coordinate gauge, a method that entails conceptual and 
computational difficulties. Here we demonstrate a gauge invariant approach to such calculations. 
For a specific set of models (requiring head on collisions and quadrupole dominance of both the first 
^\ ' and second order perturbations), we give a self contained gauge invariant formalism. Specifically, we 

^\ , give (i) wave equations and sources for first and second order gauge invariant wave functions; (ii) the 

0^ ' prescription for finding Cauchy data for those equations from initial values of the first and second 

fundamental forms on an initial hypersurface; (iii) the formula for computing the gravitational wave 
power from the evolved first and second order wave functions. 

D 

I. INTRODUCTION AND OVERVIEW 
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In the next few years gravitational wave antennas will go into operation with the possibility of detecting 

^ \ astrophysical sources. A plausible, and certainly fascinating, origin of such waves would be the powerful burst of 
■ radiation generated in the merger of two approximately equal mass black holes to form a single final hole A 
' throrough understanding of this problem will require numerical relativity on supercomputers, and is still several 
i years away. In the absence of numerical answers, some useful insights have already been supplied by "close limit" 
' perturbation theory an approximation method in which the spacetime of the merger is considered to be a 

0^ , perturbation of the spacetime of the single final hole. Close limit calculations, to first order in some separation 
parameter, have proven to give excellent agreement with numerical relativity in the case of head on collisions, models 
simple enough to be computed with numerical relativity. In principle one can apply the close limit method to collisions 
that are still beyond the scope of numerical methods, and in fact, one such result has already been reported [ p^ . 
^ \ In principle, first order perturbation theory works in the limit that some expansion parameter vanishes. In practice, 
bJO the calculations are reasonably accurate for a range of that expansion parameter, up to some maximum. A shortcoming 
^ of first order perturbation theory is that there is no indication, internal to the method, of the size of that maximum. 
• ^ , For this reason, second order close limit theory was developed jT^-p^t. When perturbations become large enough 
that the predictions of first order calculations differ significantly from the predictions of second order, it is a sign 
that perturbation theory is at its limit of applicability. Comparison with numerical relativity results, where they are 
available, has demonstrated that this method of determining "error bars" on perturbation theory is quite reliable 
[p^-JTof. Perturbation calculations must deal with the freedom to redefine coordinates, that is to do coordinate "gauge 
transformations." This can be approached in two very different ways. One way is to eliminate the coordinate freedom 
by fixing the coordinates. This, in fact, was the way in which the perturbation work on nonrotating holes p^-[l9|] 
was carried out. (Except for very recent work pO[| , these second order calculations were done as perturbations of 
a nonrotating black hole.) The even parity perturbations that were the focus of these calculations, were done with 
a second order extension of the Zerilli |^l[| formalism. In the extension to second order, the Regge- Wheeler |2^] 
conditions were used to fix the coordinates to second order as well as first. This gauge- fixed approach simplifies 
some of enormous complexity of the equations that arise, but a price has to be paid and paid twice. Explicit gauge 
transformations have to be performed to convert the initial value solution to the Regge- Wheeler gauge and an explicit 
transformation has to be performed to relate the computed perturbations to an asymptotically fiat gauge in which 
information can be extracted about radiated power. 

A second way of dealing with gauge freedom is to find combinations of metric perturbations that are gauge invariant, 
and to work only with gauge invariant quantities. The details of this procedure in first order calculations were given 
by Moncrief [ p3| , who constructed a gauge invariant combination that we shall call the Moncrief first order invariant. 
For vacuum perturbations, this Moncrief invariant satisfies the same equation as the Zerilli wavefunction and in the 
Regge- Wheeler gauge, for vacuum perturbations, the Moncrief invariant reduces to the Zerilli wavefunction. 

Moncrief 's choice is not the unique gauge invariant that can be constructed from even parity metric perturbations. 
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but it is the unique choice constructed entirely from "data" on an initial hypersurface (i.e., the 3-geometry and the 
extrinsic curvature). This property makes it especially convenient for use in calculations of evolution of perturbations. 
In work of this sort we start with a solution on an initial hypersurface in some gauge (whatever gauge is convenient for 
solving the initial value problem). In the Moncrief formulation, the starting value of the wave function to be evolved 
can be evaluated in any gauge, hence it can immediately and directly be constructed from the initial value solution. 
A gauge invariant formulation is similarly convenient for extracting information about radiation energy carried in the 
perturbations. In principle this requires that one examine the perturbations in an asymptotically flat gauge. With a 
gauge invariant formulation, the form of the invariant can be related to an asymptotic gauge in a relatively convenient 
way. 

The goal of the present paper is to cast the problem of second order perturbations of the Schwarzschild spacetime 
into a form that has the same convenience that the Moncrief formalism provides in first order perturbation calculations 
of evolution. That is, we will provide a second order formalism in which computations are carried out only with gauge 
invariant quantities, and in which these quantities are constructed from the first and second order perturbations of the 
spatial metric and of the extrinsic curvature for an initial hypersurface. The motivation behind this is to demonstrate 
the potential advantages of such a reformulation. To make this demonstration clear we choose to focus not on the 
general problem of perturbations of the Schwarzschild spacetime, but rather on a restricted class of problems for 
which a formalism can be presented with explicit details. 

One restriction is that we will consider only the axisymmetric collision of nonspinning holes. Inherent in this restric- 
tion is the simplification that odd-parity perturbations vanish at every order. The decomposition of perturbations 
into multipoles will be used, and for simplicity we will present results only for the second order quadrupole (i.e., 
£ = 2) perturbations. This choice is justified by the fact that most of the radiation is expected to appear in the 1 = 2 
multipole. 

A very different sort of restriction is that we shall consider only £ = 2 first order perturbations. In principle, first 
order multipoles of many different orders can couple to the second order quadrupole through the nonlinear mixing 
of first order multipoles. One justification for ignoring the contributions oi £ 2 first order terms is the example of 
close limit perturbations of collisions starting with conformally flat time symmetric initial data, such as the initial 
data of Brill and Lindquist |2^, or especially the solution of Misner |^^, the first and clearest example to which 
second order analysis has been applied ip^Hlq] . For these initial data sets in the case of an axisymmetric collision of 
two equal mass holes, the first order perturbations are purely quadrupolar. For other initial value solutions we might 
expect quadrupolar first order perturbations to be larger in some sense than other first order multipoles, but this 
is an inadequate justification for ignoring other first order multipoles. The real justification then is to simplify the 
presentation of very lengthy expressions that illustrate a more generally valid approach. 

In the remainder of this section we sketch out the basic ideas behind the construction of a second order gauge 
invariant; details will be given in the sections that follow. 

We consider that we have a parameterized family of spacetime metrics of the form ga^(x^, e), and that these metrics 
can be expanded as 

5«/3(x^)=5S(x^) + e5S(^") + ^£'5S(^") + --- , (1) 
where the background metric g^^^ is the Schwarzschild metric in our case; g^^ is called the first order perturbation 

(2) 

to the metric; g^^ is called the second order perturbation to the metric, and so forth. Let us now consider a 
parameterized family of coordinate transformations, also called transformations of the coordinate gauge, or simply 
"gauge transformations," x!^^^ = F'^{x°',e), that can be expanded as 

<,^=x'^-f6^W^ + ie2^(2)^ + -.. (2) 

Such a change of coordinates will transform the metric perturbations 5^^^, <?i^^, • • •. It is useful to consider special 
cases of the general transformation (|^). If ^(i)^ = 0, we call the transformation a purely second order transformation. 
Note that g^^p is invariant under this type of transformation, but gl^^^gl^^ ■ ■ ■ are not. 

The Moncrief formalism is based on a certain linear combination of first order perturbations g^^^ that can be 
determined purely from hypersurface information, i.e., from the first and second fundamental form of a hypersurface 
that, to zero order in e is a constant time surface in the Schwarzschild geometry. Moncrief shows this combination to 
be gauge invariant [invariant under transformation (^], to carry all the first order gauge invariant information, and 
to satisfy a simple wave equation, the Zerilli equation. We use to denote Moncrief's combination of first 

order perturbations g'^l . [We will present this combination explicitly below in ( |l9|) after we have introduced multipole 
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decomposition and the appropriate notation.] We will use L^"^^ to represent the same combination of second order 

perturbations g^^^- It follows immediately that the second order combination L^"^^ is invariant under purely second 

order gauge transformations. Since is constructed from the first order perturbations of the spatial geometry and 
extrinsic curvature of a hypersurface, it follows that L*-^-* can also be constructed from hypersurface information. We 
next turn to the question of the wave equation satisfied by L*-^-*. 
The vacuum Einstein equations can be written as 

E{g^,) = ^(5^ + e^S + C +■••)= (3) 

where E represents the set of nonlinear differential operators that generates the Einstein equations. The terms in 
that are zero order in e will be satisfied automatically because g'f^, the background metric, is a solution to the 
vacuum Einstein equations. To find the equations satisfied by the first order perturbation, we expand (^) in powers 
of e, and write the set of first order equations as 

O(ffS)=0. (4) 

Since the perturbations g^^^ can only appear linearly, O represents a set of linear differential operators. 

The second order part of the expansion of (H) will involve terms linear in second order perturbations and terms 
involving products of first order perturbations. These equations can be written symbolically as 

= (5) 

In this form, the products of first order terms appear on the right. If the perturbative problem is solved order by 
order, the first order problem may be considered already to have been solved, so that the right hand side of (0) 
can be considered as known. It should especially be noticed that the operator O is "zero" order. That is, O is 
precisely the same operator that appears in the first order equations . One can view (j^) as a system of differential 
equations differing from (^) only by the presence of known source terms. We know that the first order equations can 
be rearranged into a single wave equation, the Zerilli equation, which we symbolize as 

Z(*(i)) = 0. (6) 
It follows that the equations of (Q) can be rearranged to give a single wave equation 

Z{L^^^) = SMon{g^^lg^^i), (7) 

in which the differential operator Z is the Zerilli operator. The right hand side represents the set of terms quadratic 
in first order perturbations, that result from forming the Zerilli equation for the Moncrief combination. These first 
order terms can be viewed as known once the first order perturbation problem has been solved, so is to be viewed 
as a wave equation for L^^^ with a known source. 

Though the quantity L^^^ is constructed from hypersurface information and satisfies a convenient equation, it is 
not what we seek. We have seen that it is invariant under purely second order perturbations, but it is not invariant 
under more general gauge transformations. That is, L'^' will in general change under a transformation (j^) with 
^^^■^ 0. In order to construct a second order perturbation function that is gauge invariant we must add another type 
of expression to L^^-*. Let Q^^-* represent any combination of products of first order perturbations. Since the operator 
Z is linear, for any such Q'^' the quantity 

vI/(2) = l(2) + g(i) (8) 

will satisfy an equation of the form 

Z{¥'y) = Z{L^^)+Q^'^) = 5mo„(5S,5S) + ^(Q(^^) . (9) 

The added source term Z{Q'^^^) is known once the first order problem is solved, so like satisfies a Zerilli 

equation with a known source. One of the main points of this paper is to display explicit forms of Q^^^ for which 
is gauge invariant for the general gauge transformation in (g). It should be noted at the outset that Q*^^^ cannot 
be unique, and hence cannot be unique. To see this, consider ^'^^^) to be any quadratic combination of 

terms in then define 
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(10) 



and note the following: (i) ^E'^-'^' is gauge invariant under first and second order transformations, hence 'i'f^l is first and 
second order gauge invariant if is first and second order gauge invariant, (ii) Since Moncrief 's is constructed 
completely from hypersurface data, it follows that this is true also of 5"^^^ if it is true of '^^^K (iii) Like the original 

(2) 
'alt 



invariant, the alternate invariant ^/Jil satisfies a Zerilli wave equation 



with a known source. 

The remainder of this paper will be organized as follows. In Sec. || the perturbed metric tensor and wave equations, 
to first and second order, are introduced. The details of first and second order gauge transformations are discussed 
in Sec. [II and a wave function is presented that is invariant under these transformations. The procedure for finding 
Cauchy data for this wave function is given in Sec. [V. The relationship of the invariant wavefunction to gravitational 
wave energy is analyzed in Sec, and a summary and discussion are given in Sec. VI. Throughout the paper we use 
the conventions of Misner et al. [26|. In particular we use a metric with sign conventions -H — I — h, and units in which 
c = G = 1. 



II. PERTURBED METRIC TENSOR AND MONCRIEF WAVE EQUATION FOR SECOND ORDER 

A. Perturbed metric tensor 



As was discussed in Sec.[y, to both first and second order we consider only i — 2 multipoles. Using the standard 
Regge- Wheeler ||2^ notation we can write the perturbations [the g^^poi (|^)] to the Schwarzschild metric for mass Af , 
as 



gtt = -(1 - 2Af/r) 

gtr 
gte 

grr = (1 - 2M/r)-^ 



eH['\r,t) + '^H['\r,t) 
S^\r,t) + '\'i\r,t) 



P2{0) 



l + {eH^'\r,t) + '^H^^\r,t)\P^{e) 



(12) 
(13) 
(14) 
(15) 
(16) 
(17) 
(18) 



Here we are using P2{0) to denote the Legendre polynomial of order 2, with argument cos 6*. By P2{d) and (^) 
we mean respectively the first and second derivative of P2(0) with respect to 0. We include an upper index between 
parentheses whenever it is necessary to clarify the order of a perturbation quantity. 



gro 



eh['\r,t) + i-h^^Hr,t) 



+ leK^^Hr, t) + ^ii:(2)(r,t)lsin2 eP2{e) + {eG''^\r, t) + ^^G^'^\r,t) \ sin 61 cos 6IP2 (6*) 



B. Zerilli wave equation 

As described in Sec. |, a starting point in the search for second order invariants is the first order Moncrief invariant 
For £ = 2, in terms of the Regge- Wheeler notation introduced above, is the following linear combination 
of first order perturbations 
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6 (2r + 3M) 



2{r - 2M){H^ ' - 



(19) 



Moncrief has shown that this first order combination is invariant under gauge transformation (||) and satisfies the 
ZeriUi 111 equation 



Here the Zerilli operator Z is 



Z{¥^^) = . (20) 



where r, is the usual "tortoise" coordinate covering the exterior of the black hole, 

r, =r + 2Mln(r/2A/-l) , (22) 

so that the horizon is at = — cx) and spatial infinity is at = cxd. The potential term in the £ = 2 Zerilli operator 
is given by 

^('^^ r3(2r + 3M)^ ' ^^'^ 



We now define L*^^^ to be the second order equivalent of L^^^ 



6 (2r + 3Af ) 



2(r - 2M){Bf^ - TdrK^^y) - 2{r - 3M)K^^^ + 6 \rK^^^ + il-^Ml (^23^^(2) _ 



(24) 



As was argued in connection with (^ , this second order combination satisfies an equation of the form (|^) with S'Mon 
a sum of products of the first order perturbations, Hq^\ h['^\ h'^\ h'^\ h^i \ K^^\ G^^\ and the derivatives of these 

functions. The explicit form of Syion is straightforward to compute; one simply repeats the steps that lead to the first 
order Zerilli equation (^0|) and keeps all terms of second order. But the result is extremely lengthy and will not be 
displayed here. 

It is worth noting that Syion is automatically invariant for purely second order transformations since first order 
perturbations do not change for purely second order transformations. Since Z is invariant under general gauge 
transformations, (j^) then tells us that L^^-* must be gauge invariant under purely second order transformations. This 
is a property that also follows from the manner in which L^^^ is constructed. The validity of (Q) can then be viewed 
as a check of consistency. 

III. SECOND ORDER INVARIANT WAVEFUNCTION 

To explore the gauge changes in metric perturbations we must introduce a specific form of a gauge transformation. 
The form of the transformation in (^) was used in the studies by Gleiser et al. [p^-|l6| . Here we choose instead the 
equivalent form of higher order gauge transformation given by Bruni et al. p7[ | 

x<^^^ = + e e(i)"(a;^) + {e^ 12) [^('^"(x'^) + a^C^'^^l^") ^^'^''(x'')] ■ (25) 

The reason for this choice is a practical one. With this notation ^(i)" and ^(2)q ^^^^ treated as generating vectors 
and the first and second order gauge transformations can be written in the form of Lie derivatives 

(26) 

= (/:^(., + /:|<,) ) ffS + 2/:5(i, • (27) 

Lie derivatives of tensors can be handled automatically by the Maple symbolic manipulation language that was used 
to do the computations. 
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Since we are using multipole decomposition and keeping only the quadrupole terms, we write the components of 
the generating vectors as 



The first order gauge transformations then take the form 



Sh^n 



2drCi^\r, t) - (1 _ 2M/ry^ 2MCi^\r,t) 
(f - 2M/r)-^ Cl^\r, t) + r''drd/\r,t) 

2C«(r,t) 



(1 - 2M/ry^ dtCi^\r,t) - (1 - 2M/r) drCy{r,t) 



{l-2M/r)cl'\r,t)+r'dtCi'\r,t) 



(28) 
(29) 

(30) 
(31) 
(32) 
(33) 
(34) 
(35) 



Pure second order gauge transformations would look exactly like (pOD-(p5D, replacing upper index 1 by 2. An example 
of general second order gauge transformations, after projecting into ^ = 2, is 

JG(^) = 2Cf ) + ,^,l'\,^ i^r^ir ^ 2M)C^^d,&^^ - 

3r2c(i)C^i) - 18r3(r - 2M)&^^d^^ + 2r3(r - 2M)Ci^'^ drG^^'^ - 

6r(r - 2M)/i^^'c(i) - 6r(r - 2M)h\^^&^'^ + 3(r - 2Mfc[^^c[^'^ + 2r3(r - 2M)c[^^ dtC'^^^ + 
Ar^ir - 2M)&pK'^^^ - 18r3(r - 2M)C'^^^ G^^'^ + 

4r2(r - 2M)C(i)g(i) -^[r- 2M)Ci^^ drC^''^ + Sr^r - 2Myf^C'^^^) . (36) 

Expressions like (pq) generate the first order gauge transformations of the linear part in second order perturbations 
ofL(2). ^ 

The key to building invariants at second order is to take combinations of the gauge transformation equations in 
(p0|)-(p5|) that isolate the coefficient functions occurring in (28) and (p9). 



(r-2M) 



(37) 

(38) 
(39) 



With these at hand one can construct first order quadratic terms that, under a first order transformation, cancel 
gauge dependent terms arising from the transformation of L^'^\ This procedure leads to the following as the simplest 
choice for a first and second order gauge invariant. 



*(2) = r (2) , Q 
RW — ^ W 



(1) 
RW 



(40) 



where L^^' is the second order equivalent of 'i>^^\ given in (p3), and 



{r - 2M)2(/iW _ !_9^,G-(i))(_2r(3r + M)K^^^ + 6r(r - 2M)drh[^h 



21r3(r~2Af)(2r + 3Af) 
6r(5r + 4M)G(^) - Ir'^MdrK^-^^ + r(3r + 2M)H^^^ - 3r^(r - 2A//)a^2G(^) - 6(5r + 4M)/i^^' + 
Qr^{r + bM)drG^^^ + ^{r - 2M)dr2K^^^ -r^{r~ 2M)drH^^'^) + 

r\r - 2M fd/-^—^ih['^ - ^9.G(i'))(-2rijf ^ + 6h['^ - Sr^drG^'^ - 6rG^'^ + v^drK^^^ + 2rif (D) + 



r'^ir ^ ^^'^) g(i)(i8(y _ 2M)h^l^ - 2,r'^{r - 2M)drG^^^ - 3r^{r - 2M)drK^^'^ - 24r(2r + 3M)G(i) + 
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9r(2r + 3Af)i^(i) - 3r(r - 2M)Hi'^) + _ _ 2M)h^^^)drG^'^ + 

3r3(r - 2M)artG(i) + 3r(r - 2A.f)af/i^^^ + 3r(r - 2M)h[^'> + 3r(r - 2M)drh'^„^^ + 6r^{2r + ZM)dtG^^^ - 
6(5r + 4M)4^^) + 

r2(r - 2M)2a,(^-izl_(4i) _ !^a^G'(i)))(-3r2a,G(i) + G^'^ - 2(r - 2M)i7W + r^dtK^^^) + 
8(r - 2M)V^(/^« - !^a.G(i))a.(^^-^(/.« - ^9.G(i))) - 

2 

6r(r - 2M)2(5r + AM){h^^^ - y c'rG(^))G(^) + 



3r2(r - 2M)3(/i(i' - a^G(i))a.G(i) 



r%r-2MYdr{ 



2^ ,(r-2Af)^ .1) 



-a.GW))9.(^^l^^(/^« - ;^a.G«)) 



6r^(2r + 3M)(r - 2M)G^^'^G'-^^ - 3r^{r - 2MfdrG'^^^G'-^^ + ^r^{r - 2M)3a^G(i)arG(i) 
3(5. + 3M)(. - 2M)^(^-l^(/.« - ... _ ^^^^ 



5*G«))(r^47T(4''-^a,G«)) 



(41) 



In the Regge- Wheeler choice ol gauge the perturbation functions h'^\ /i^^-* and G^^^ are chosen to vanish. In this 
gauge the right hand side of ( pj ) vanishes. For this reason we denote by the particular choice of quadratic terms 
appearing in (|4|) and we denote the corresponding wave function as ^'^w- Note that ^'rw becomes simply L^^^ in 



(42) 



(2) 

the Regge- Wheeler gauge. The wave function ^^i^ satisfies a wave equation 

Zi^'il,) = Z{L^^^ + g|^^) = 5Mo„(ffS,5S) + Z{Q^^l,) ^ 5rw 

If one works explicitly in the Regge- Wheeler gauge, then the terms in Z(Qp^) vanish. 

The expUcit expression for the source term is extremely lengthy. To compress it into a manageable form we introduce 
a few simplifications in notation. We use a prime (') to denote partial differentiation with respect to r and a dot (') 
to denote partial differentiation with respect to time, and we write it in terms of two first order quantities K and H2'. 



'RW 



-2/189) (1512M^ - AOQ^r^M"^ + 1602r''Af - 252.^ - 4248rAf'^ + 5A90r^M^)H2K' 

2Mr^ - %72r^M + 324Af* - 468rA//3 + llQlr'^M'^)KK + 

30r^ - ISQr^Af* + Mlr^M^ - 77r^M^ + 432rM^ - 52r^M)kk + 

-9r^ + 72r^M - 216r^Af2 + 288r3Af3 - 144r2Af^)i^'if' + 

-36r^ - 864r^Af2 + 288r^Af + llb2r^M^ - 576r^M'^)KK" + 

729r^Af2 - 459r^A// + 72r*^ + 882r^M^ + 2376rAf ^ - 3348r2Af^)i?2-ft:' + 

-783r''Af2 + 90r^A/ + 2538r^A/^ + 1944rAf^ - 3636r^M*)H2K" + 

-306r*M + 1944r^Af 2 - 864A//^ + 3744rAf^ - 'i320r^M^)H2H2 + 
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iSTSr^M^ + 387r^ - 1359Mr^ + 3132M^ - 2376rM3)i?2-f?2 + 

'-4:50r^M + 1080M'' + 2070r^M^ - 2880rM^)H2K + 

-36r^ - 87/m2 + ASr^M + 5'i0r^M^)H2K + 

36r^ - 276r*M2 + 162r^M^ + 540r'^M'^ - A2r^M)Kk + 

-lllr^M + 18r^ - 576r^M^ + 216r^M^ + 72r^M^ + 432r^M^)A''i?2 + 

-mr^M + 18r^ + 72r^M^ - 576r^M^ + 432r^M^ + 216r^M^)H2K"' + 

'-56r^M^ - 2r^ + ASr'^M'^ - 3r^M - 96r^M^ + UAr^M^)kK"' + 

-bAr^M + 12r^ + 12Qr^M^ - lUr^M'^ + 36r^M'^)KH2 + 

-36r^M'^ + hAr'M - 12r^ + lUr^M^ - 120r^M^)Kk" + 

-eOr^M^ - 18r^M2 - 6r^ + 27r^M + 72r^M^)K'k" + 

-ISr'^M^ - 6r^ + 27r^M + 72r'^M'^ - 60r^M^)K"k' + 

-48r^Af2 - 16r^M^ + 2Qr^M - Ar^ - mr'^M^ + I28r^ M'^)k' H'^ + 

iGr^M^ + 4r"' - 128/' Af^ - 26r^M + 48r*Af2 + 9Qr^M^)k'k"' + 

[Sr'^M - ASr^M^ + 2r^ + QQr^M'^ - 144r-^A'/-^ + b&r'^M^)kH'^ + 

IhlrH'I - 258r^M^ - 22r^ + 1288r^M'^ - 960r^M^ - 316r^ M^)k' H'^ + 

;i8r^ + 588r^M'^ - 720r^M^ - 133r^M^ + 12Ar^M^ - 29r^M)H'^k + 

-108r^ + 792r''Af - 2214r3Af2 + A32M^ - 1800rM'^ + 2'dlQr'^M^)KK' + 

l7lr^M - 24r^ + 2376r2Af^ - 1362r^Af^ - 211r^M^ - 97r^M'^)k'k + 

'57r^M - 6r^ + 378r^M^ - 99r^M^ - 135r^M'^)K'k + 

A32M^ - 36r^M + 162r^M'^ - 360rM^ - 108r^M^)H'2K + 

'-18r^ - 72r^M^ + 117r^M - 216r^M^ - 432rM^ + 576r'^M^)H2H2 + 

36r^ + lAAr^M^ - 23Ar^M + A32r^M'^ - 1152r^M* + 86ArM^)KH2 + 

'-128r^M'^ + 4ri° + Wr'^M^ + A8r^M'^ - 26r^M + 96r^M^)k"k" + 

leOr'^M^ - 17r^ - Gr'^M + 1992r^M'^ - 226Ar'^M^ + A7Ar^M^)k'k' + 

'Ar^ - 26r^Af + 9Qr^M^ - l28r^M^ + A8r^M'^ + IQr^ M^)H'^H'2 + 

'A8r^M'^ - 2Qr'^M + Ar^ + Oer^A/^ - 128r''Af'' + l&r^M^)H2H'^ + 

-IQr^M^ - 48r^Af2 - 4r^ + 2Qr^M ~ 9Qr'^M^ + l28r'^M'^)H2k"' + 

-18r^Af2 + 27r^Af - 6r^ + 72r'*Af'' - GOr'^ M^)H'2k' + 

[-bAr^M + 12r^ - 144r3A'/'' + 120r*M^ + 36r'^ M^)H2H'2 + 

'-57r^M + 6r^ - 252r^Af^ + UOr^M^ + 93r^M^)K' H2 + 

-3Qr^M'^ + UAr^M^)H2k" + 

Ulr'^M - 18r^ - 816r^M^ + 1120r*Af* - 192r^M'^ - 312r^M^)H2k" + 
'33r'^M - 6r^ + lOSr^Af* - 72r^Af^ - 33r^M^)K"k + 
ISOr^M - 24r^ - 384r*A4'^ + 720r^Af^ - 252r^M^)Kk' + 
'-27r^M + 6r^ - 72r^Af^ + 60r'^M^ + I8r^ M'^)h'2H2 + 

-144r^Af2 + 117r^A/ - 18r^ - 252r'^Af^ + A32r'^M'^)K' k' + 
-32r^Af^ - 96r^Af2 + 52r^Af - 8r'^ - 192r*M^ + 256r'^M*)H'^k" + 
232r^M^ + 384r^Af2 - 179r*^Af + 22r'' + 1296r'*A^^ - 1568r'^ M*)k' k" + 
182r^M^ + 43r^Af - 18r^ + 1224r-''A/^ - 672r'^Ai"' - 362r^M^)kk" + 
Mr'^M - 12r^ - l2Qr^M^)H2k" + 
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-21r'^M + 6r® - 72r'^M'^ + eOr^M^ + \^r^M^)K" H2 + 

60r^ + 2604r3Af'' - 1584r2M5 + GTr'^M^ - QSSr-^M^ - l?,?,r^M)H2K' + 

12r^ + lOSr^M"* + iQr^M^ ~ 69r^M^ - l2r^M)H2k + 

-42r^ + 720r3M'' - 492r^M^ - 21Qr^M'^ + 225r^M)H2K' + 

-24r^ - TSr^M^ - 120r^M + 168r*M^ - 576rM^ + lHQr"^ M^)H2K + 

36r^ - llb2r^M^ + 144r''M^ + Ai2r^M^ - 234r^M + 8Mr^M^)H'2K" + 

-36r^ - 144r3M3 - 432r''Af2 + llb2r'^M^ - 8MrM^ + 234r^M)i/2if2 + 

ISr'^M^ - 27rHl + 6r^ - 72r^A/'' + mr^M^)K' H'^ + 

18r^ + 480r2M^ - 840r3M^ + 396r''M3 - 93r^M + mr^M'^)H2H'^ + 

iUr^M^ + 42r^M2 - 195r^M + 54r^ - 360r'^M*)H2H2 + 



-eeOr^M* + 364r-^Ar - Qlr'^M^ + 94r^M - 33r'' + 192rM^)ff2H2 + 
150r^M^ + 120r^M^ - QOr^M - 360r^M^)KH2 /r^(2r + 3Mf{r - 2M) . 
The quantities K and _ff2 occurring in (Bsf) are the foUowing combinations of first order metric perturbations: 



K = K 



+ (r - 2M) {drG^^^ 

+ (2r - 3M) (^5^G(i) - ^/^l^^) + r{r - 2M)dr {drG^^^ 



(43) 

(44) 
(45) 



Note that K and H2 reduce to the metric perturbations K^"^^ and in the RW gauge. 

As explained in connection with (|lO|), it is easy to construct an alternative wave function by adding to terms 
quadratic in Any such modification, however, will not have the property that quadratic terms vanish in the 

Regge- Wheeler gauge, and might make the expression for the source more complicated. 



IV. INITIAL DATA 



A. Unit normal to the perturbed hypersurfaces 

Our spacetime will be foliated, outside the horizon, by spacelike surfaces of constant coordinate time t that, to 
zero order in e, agree with the surfaces of constant Schwarzschild coordinate time. Our initial hypersurface will be 
assumed to be one of these constant t surfaces. Viewed as a foliation of a given spacetime, the constant t surfaces will 
change under coordinate transformation, but our variables ^^^^ and ^r^j specified at a particular coordinate value, 
are invariant under such transformation. More subtle is the fact that for a fixed foliation of a fixed spacetime, the 
meaning of the partial derivative d/dt will change when gauge transformations are made that perturbatively change 
the spatial coordinate labels on a hypersurface. But our quantities are invariant with respect to such perturbative 
diffeomorphisms on the hypersurfaces, so the d/dt operation on our quantities will be unchanged. 

To relate the time derivatives of our quantities to the extrinsic curvature we must first find, to first and second 
order, the components of the future directed unit normal to the constant t hypersurfaces. The computation starts 
with the relationship of the metric and the shift components Ni in the standard ADM p3| decomposition 



(46) 
(47) 
(48) 

With the definition of the perturbed metric tensor (|12|)-(|18|) and the definitions for the shift (^)-(^8|), the perturbative 
expressions for the shift components given to second order are 



gtr = Nr{r,t) 
gte = Ng{r,t) 
9t4, = N^{r,t) 



Nr{r,t)= eH['^ + ie'/2)H^'> P2{e) 



(2) 



Ne{r,t) = 



N^{r,t)=0 . 



,(1) 



(6V2)4'^ P2{9) 



(49) 

(50) 
(51) 
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The vanishing of the shift component A''^ is due to the assumed axisymmetric character of the colhsion. 

In order to obtain the expression for the perturbed lapse to second order we can use the ADM lES] equahty for 

(4), 



'9tt, 



(52) 



Here summation on i and j ranges over r, 9 and (f>, and ^^"f g'^^ is the inverse of the 3-metric on a constant t hypersurface. 
Using (^9|)-(pll) in (|5^) gives the perturbed lapse to second order 



N{r,t) = - 2Af/r 



hi'hl'^P,{9)P,i9) , 1^(1)^(1) 



r{r~ 2M) 



H^'>H^'>P,{9)P2{9) 



(53) 



In a similar manner one can compute and to second order from iV* — ^'^^g^^Nj. From these, and TV, the 
covariant 11^ and contravariant n" components of the orthonormal vector to the time slices are given as: 



n„ = (-7V, 0,0,0) 



(54) 



1 

n° = (—,- — ,- — ,0) 



(55) 



B. Multipole decomposition of extrinsic curvature 



In computing the extrinsic curvature Kab on our constant t hypersurfaces, it is useful to start with a tensor harmonic 
decomposition, similar to that for the metric tensor given in (|l2|)-(|l|). We define the quantities xi]} , kI^J , and 
K^"^ by the relations 



Krr 

Kee 
Km 



(1 - 2M/r)- 



^-K^K^r) + -K\^>{r) } P,{9) + \ eK^^>{r) + -K\^>{r) \ P, (9) 



'(2), 



P2i0) 



'(2)/ 



^K^xHr) + -K'^>{r) \ sin' 9P^{9) + <j eK'^>{r) + -K)^\r) \ sm9 cos 9P^{9) 



'(2), 



(56) 
(57) 
(58) 
(59) 



C. Time derivatives of the metric perturbations 



To relate the extrinsic curvature and time derivative of metric perturbations we will need the projector 

/f = (5f + n.n'' . (60) 
into the 3-dimensional hypersurfaces. The extrinsic curvature is defined, in terms of the unit normal n", as 

K^, = -IPtP^n^f.;.) ■ (61) 

Here the symbol "; " means covariant derivative with respect to the 4-metric tensor. With the relations in (p6|)-(p5|) the 
right hand side of (^l|) can be decomposed into multipoles with coefficients expressed in terms of metric perturbations. 
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With ( pl\ ) these coefficients are then related to the multipole coefficients defined in (^6|) - (59). As examples of the 
results of this procedure, the time derivatives of TJj"'^-' and iJj^-* are 



dtH^^^^ = 2(1 - 2M/r) drH[^^ + {2M/r^) h[^'> + 2^1- 2M/r K^^ , 



(62) 



12Mh\''ho 



■tn^ = 2(1 - 2M/r) drH[^' + (2M/r2) ijf ' + 2^1 - 2M/r K^^^ - {2/7r^) 

12r^hi''>drh[''^ ~ 2Arh^^^drh^l^ - 2r^ H^^^ OrH^"^ + 2r^ H^"^ drH^^^ - ArHl H^^^ drH^^ 



r(l) 



„4 rr(l) 



(1) 







(63) 



These examples show the general pattern. Time derivatives of first order perturbations of the metric tensor are 
expressed as linear combinatins of first order metric perturbations. The expressions for time derivatives of second 
order perturbations contain terms linear in second order metric perturbations, and terms quadratic in first order 
perturbations. 



D. Extrinsic curvature and momentum constraints 

Since axisymmetry is assumed, the nontrivial first order momentum constraint equations are R^l) = and R^g^ ~ 0. 
The partial derivatives, with respect to time, of the first order metric perturbations can be reexpressed in terms of 
components of extrinsic curvature by using relations like (^2|). The resulting, simplified, momentum constraints are 

= r^{r- 2M) kI]} -r'^{r - 2M) K^^^ +r^{r- 2M) K^^^ + {-2r'^ + 7rM - 61/^) -r^{r - 2Mf drKlg' (64) 
^ -r^{r - 2Mf drX'^^^ + ir^{r - 2Mf drK^^^ + r{r - 2Mf K^l^ - r(r - 2Mf k''^^ + 3r(r - 2Mf k''^' 

-3(r-2Af)2ir« . (65) 

These results will be needed below to help connect dt^^m ^"^^ initial value data. 



E. Wave function on the initial hypersurface 

Cauchy data for the first and second order wave equations requires the values of 4"*^^^ and on the initial 

hypersurface. If the 3-geometry of the initial hypersurface is known, ^'^^^ follows immediately from the definition in 
(p^), and L^^) similarly follows immediately from (p4|). But the specification of ^'^t^^ requires Q^^' ^"^"^ ^'^^ definition 

in ( ^l|) appears to require a knowledge of terms that do not follow directly from the hypersurface 3-geometry. 
All such terms can be grouped into only two combinations of perturbations: 

[4'^ - Y5,G(i)] and [H^"^ . (66) 

These two first order expressions, however, appear in the following components of the extrinsic curvature: 

^^^^ = [ (4'^ - t5,G(i))P;((?) + (r - 2M) - '\.jM ^'^)P2{9)] (67) 



= ^===[sin0cos0 (/.W - ^a,G(i))P;(0) + (r - 2M) sin^ (i/f ' - -^J^^dtK^^))P^{e)] . (68) 



An initial value solution in Einstein's theory consists of both the 3-geometry and extrinsic curvature of the initial 
hypersurface. From an initial value solution, then, the terms in ( |66| ) can be evaluated, and the process of specifying 

^''^^^ and 4*^1^ '^^^ t)'^ completed. 



11 



F. Time derivative of the wave function on the initial hypersurface 



The complete specification of Cauchy data for tlie first and second order Zerilli equations includes the time deriva- 
tives 5t*(i) and at*[^^. These require the time derivatives of metric perturbations, which are found starting with 
(|6T|), as shown in examples ( p2[ ) and (|6^). With this approach the computation of St^'^^-' is straightforward. The 
computation of St^'j^^, however, requires dfQ^^ and is not straightforward. The evaluation of this time derivative 

produces terms that involve H^^^ multiplied by time derivatives of groups of extrinsic curvature terms. These groups 
of terms turn out to be those that occur in the momentum constraint (|6J), so that the troublesome term is guaranteed 

to vanish. The resulting, simplified expressions for St^''-^^ and 9t^'p^, in terms of hypersurface information, are 



(1) 



21r4(2r + 3M)y^l ~ 2M/r 
(28r'*Af2 - 28r''^M + 7r^)K^^l^ + (168r''M - 168r^M^ - 42r^)i^^g^ + 



(69) 



(2) 



t*RW 



(28r*'M - 28r'W^ - 7r')drK'^^ + {-SAr^M + 2lr' + Mr''NP)drK)^ 

(2) 



(2) 



21r4(2r + Md)^! - 2M/r 

28r^M2 - 28r^M + 7r^)K^f + (168r''Af - 168r^M'^ - 42r^)X^g^ + 



14r^ - 28r^M^ 
4r^ - 32r^M - l28r^M^ 
4r^ - 32r^A//3 - 24r^A'/ - 



'(2) 



-2r^ + 12r^M 
8r^ + 216r^M2 



16r^A/3 - 24r^M2)a^2G(i)9rXjf^^ + 



14r^ + 22r^M + AGAr'^M^ - 2Mr^M^ - \h2r'^M^)d^2G^'^^ 
-Ar^ - 96r^Af2 + 32r^Af - 64r^Af'' + \28r^M^)dr2G'^^^ drK, 



(1) 



(1) 



(1) 



-2r^ - 24r^A/2 + 12r*A/ + 16r^M^)drK'^^^dr2K^ 
2r^° - 64r^A/3 + 32r^M'^ - IGr^M + A8r^ M^)dr2G^^^ dr2Kl^^ + 
-2r^ - Ar'^M - A8r^M^ + AQr^ M'^)dr2K^^'^ k[^^ 



G 

-Ar^ + 32r^A/ - QAr^M^ + 128r^Af3 - mr^ M'^)drh'\^^ dr2K^(^ 



,(1); 



'(1) 



-2r^ ~ 24r^A/2 + 12r*A/ + l&r^ M^)drG^^^ dr2K'^^^ + 
2r^ ~ IQr^M^ - 12r*Af + 2Ar'^ M^)dr^G^^^ K^^^ + 
-2ri° + 64r^Af^ + 16r^Af - 48r^Af2 - ?,2r^ M'^)dr<iG'^^^ drX''^^ + 
-112r*^Af^ + 120r^A/2 _ 52r**Af + 32r-'^Af'' + 8r^)drG'^^^ dr2K''^^ + 
6r^Af - Sr^Af^ _ 2r^)dr2G^^^ k''^^ + 
-24r^Af - 32r^Af''^ + 48r'^Af2 + Ar^)dr2G^^^ K^^^ + 
2Ar'^M + ■i2r^M^ - A8r^M'^ 
~2ArHl + 48r''A'f 2 - 32r^Af '"^ 



Ar'^)H^^^dr2K, 



(1) 

G 



A8r^M^ - mr'^M^ - 2Ar^M + Ar'^ )h^^'^ 8^2 K 



'(1) 

K 
,3 A,ri\ 



-48r5A/2 - 8r^ - 32r*A/3 + 40r^Af + QAr^ M^)h^^^ dr2K''(^' 
-8Qr'^M^ + 96r3A/^ + 36r^A/ - 8r^ - 2Ar^ M'^)dr2G^^'^ k''^J + 
2Ar^M - A8r^K'P + 32r^Af3 - Ar'^)dr2h''^^ K^^^ + 



'(1) 
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96r^M^ - 48r''M2 + 8r^M - 6Ar'^M*)dr2G'-^^drK^G^ + 

'56r^M^ + 30r^M - 4r^ - 72r^M^)drG^'^^drK^j^^ + 

'2Ar^M - 6r'' - 2Ar^ M'^)G^^^ drK^^'' + 

-lAAr'^M^ - 14r^ + lOAr^M^ + 12r^M)drK^^'> K^^^ + 

-lUr'^M ~ 32r^M^ + UOr'^M^ + 96r^M* + 40r^)/if ^a^iiT^^^ + 

-IGr^M + 64:r^M^ - 64:r^M^)drK'^^^ K^.^^ + 

'-24r'^M^ + 2Ar^M - 6r^)G^^^ K^/^^ + 

'2Ar'^M^ - 28r^M'^ ~ 4r'' + I6r^ M)drG'^^'> K'j^'> + 

dO'ir'^M^ + 58r^ - 64/' Af - 356r'^ M'^)drG'-^'> K^^^ + 

;96/Af2 - 48/M - 64:r^M^ + 8/')K^^^ drK^^J + (SOO/M^ + 42/M - mr^)G^^'> K^^^ + 

12/ + 96/M2 - 48/M3 - 60/Af)/4^^ A',^.^' + 

;i6/M^ - 24/M^ + 12r^M - 2r^)drH^^^drK^a^ + 

;i36/M + 112/M3 - 24/ - 232r^M'^)h']^^ k'^j^^ + 

-ge/M^ - 8r^ + GAr^M^ + A8r^M)H^'^^drK!;.l^ + 

;192/m2 - 64/M - 256r^M^ + 128r^M'^ + 8r^)drh^l^ drK';^^ + 

;8/M^ + 8r^M'^ + 4/ - lAr^ M)drK^^^ drK''^^ + 

-2/ + 8/M^ - 16/M^ + 4/M)if^^^arAr^^^ + 

-16/M^ - 12/M + 24/M2 + 2r^)drG'^^'' drK^)) + 

;16/M + 256r^M^ + 64r^M* - 240/M^ + lQr^)drG'^^'i k''^^ + 

-144/M + 36/ + lAAr^M'^)G^^^Kf^ + 

;16/M^ + 2Qr^M - 4r^ - 32r^ M'^)K'^^'> drK'^^'^ + 

-72/M + 16r^ + 48r^M^ + 160/M^ - 192rM^)drh^P K^^^ + 

-976/ Af3 + 928/A/"2 + 12/ + 352rM* - 288r'^ M)h''^^ k'}^^ + 

Qir^M^ - 80r*M + 16r^ + 96r^M^ - 128rM*)hY^ drK^^ + 

;64/A/'3 - 24/ - 160/M2 + 112r'^M)K^^'> K^.l'' + (12/ + 112/M2 - 64/ Af^ _ 64r'^M)H^^'> K^^ + 

-88/ + 256/ AT + 96/ Af^ - 640rA4^ + 256A4*)hY^ K^^^ + 

44/ + 18/ A/ - 212/A/)if (i)ifW (_g^7 ^ 24/M - 24/A/'2)G(i)5rii:^^^ + 

36/ - 24/ AT - 96r^M^)G'^^^K'j^^ + 

-32/Af3 + 4r^ + 48/Af2 - 2Ar^M)drh[^^drK''j^^ + 

-64/ Af^ - 240/ Af2 + 224/ AT^ - 16r^ + 104:r^M)drG^^'^drKl.g' + 

288/ AT^ - 276/ AT^ - 336/ A/ + 228/Af - 60/)ar/i^^^A:^^^ + 

-32/Af3 + 48/Af2 + 4r^ - 2Ar^ M)drK^^^ OvK''^^ + 

;88/Af3 - 96/Af2 + 30/ Af - 2r^)drH^'^^ K^^'> + 

-48/Af2 + 24/ Af - 4/ + 32/Ar^)9^i?f 'if^g^ + 

-40/Af3 + 48/Af2 - 18/ A/ + 2r^)(9rG(^)ifW + 

192/ Af^ + 20/ - 108/ Af - 112/Ar3)9r/i[^^ AT^^^^ + 

;48/Af2 - 12r^Af - 48r^M^)h^^^drK^^^ + 

;28/ - 12/ Af - 88/Af=^)ijf 'if^^^ + 
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(-ITGr^M^ + 132r^M2 + 16r^M* ~ 8r^M - 8r^)drG^^^ drK^^^ + 



K 



8r^M + 2r'^ + 8r^M'^)drH^2^K^^A . (70) 



V. ENERGY RADIATED 

Radiation is most clearly analyzed in a reference system that is asymptotically flat (AF). By this we shall mean 
a coordinate system in which the deviations of the metric 5g^i, from Minkowski form decrease with r, at constant 
t — r* , according to: 

Sgtt, Sgtr, Sgrr ~ 0{r~^) 6gtr, Sgre ~ 0{r^^) 

Sgee + Sge^,/ sin^ 9 - ©(r^) Sgg4„ Sgeg - Sgg^/ sin^ 9 - 0{r^) . (71) 
In such a coordinate system, the power carried by gravitational waves is given by p3,B0[ 



dPower 



dn 167rr2 



sm 



1 f dge 



^9\dt 4\ dt sin^ 9 dt 



(72) 



In the statement of the conditions for an AF gauge, and in ([7 2D, there is no reference to the order of the metric 
perturbations in some expansion parameter e. To compute the power to first order in e, one finds the perturbations 
5gee,5ge^,5g^^ in a gauge that satisfies the AF conditions to first order in e and uses those values in (|7^). For a 
computation correct to second order, the values of Sggg, Sgg^, Sg^^ used in ( [72| ) must in principle be computed in a 
gauge that satisfies the AF conditions to second order in e. 

In practice, the use of gauge invariant quantities 5"^^^ and ^'p^ in our computations simplifies the evaluation of 
energy. A first and second order gauge transformation can always be done to bring the metric perturbations into 
a form that is AF to first and second order, and and 5*^^ unaffected by such a transformation. We can, 
therefore, treat them as if they had been computed in a gauge that is first and second order AF, and we need only 
read off the first and second order metric perturbations from and ^'pX- When restricted to first order only, this 
is the method that was used to compute radiated power in Refs. and po|| . 

If we take the expressions in (|l^) and ( p^ ) to be in the AF gauge, then for our axisymmetric quadrupole example 
the expression in ( [72| ) becomes 

The first order part is easily dealt with, and the treatment is identical to that in equations (III-21) - (III-26) of Ref. 
pO| . From the definition of ^f^^^^ in ( p9| ) and the conditions in (^Tj) it follows that in the AF gauge 

12rG(i) = +C>(r-i) . (74) 



Evaluating the second order part is considerably more difficult. Using the AF conditions in (|71|), one must solve 
( p^ ) for G^'^\ The result will contain terms that are linear in second order perturbations and quadratic in first order 
perturbations. The linear terms are identical to those in the first order analysis, so the result can be written in the 
form 



12rG(2)^*(,^-gW +0(r-i). (75) 
ier perturba 

(^). It remains to find ^'^ 0{r^) in the AF gauge. The result must be expressible entirely in terms of 



Here is the set of terms quadratic in first order perturbations, as defined in (0), and as explicitly exhibited in 
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since (^3|) is an expression for a physical quantity, and '^^^^ contains all the gauge invariant first order information 
about perturbations. 

To find in the AF gauge, we start by writing the following asymptotic expansions for first order perturbative 
quantities: 



rr(l) 


J- / 3 


H 






- h/r' 


H 




^2 


= f2/r' 


H 






= fs/r 4 


' lo/r^ H 






- h/r -i 










- /e/r^ - 








- /o/r' - 





(76) 

where /i • • • /lo, /o, -fi are functions of < — r*. When ( [76| ) is used in the expression for given in (pl|), the result is 
an expression of the form 

Q^^l^Qo + rQi+r^Q2 + 0{r-^) , (77) 

where Qo, Qi, Q2 are functions of t — r*. As an example, we present here the explicit expression for Q2: 

Q2 = -2/^/^' + 10/^/^' , (78) 

where a prime (') denotes differentiation with respect to f — r*. The expressions for Qo and Qi are very lengthy and 
will not be explicitly exhibited. 

The radiated power is a physical quantity. Since ^E"*^^^ and ^'j^ij^ carry all the (first and second order) gauge invariant 
information about the spacetime, it must be possible to express the coefficients Qo, Qi, Q2 in terms of To find 

such expressions we start by using the AF expansions of ( [76[ ) in Einstein's vacuum equations to derive 

3/8 - /5 = (79) 

6M/^ - 3/{o + 7M/5 - 2/6 - 3J( + /^ - OM/g - 3/3 -3/^ = (80) 

/5 - /8 - /g = (81) 

/3 + /4 + /g = (82) 

4/^0 - 4M/1 + 2I[ + 2/9 - 2M/8 - 4Mf^ + 2/4 + 2/^=0 (83) 

/e' - 3/9 = (84) 

2M/^ - /^ + 2/^' - 9Af /^ . (85) 



The first two equations are respectively the and parts of the Einstein equation Gu + Gtr = 0; Eqs. (|8l]), ( |8^) , 
( p3| ) are combinations of the and parts of Gee = and G^^ = 0; Eq. ( p^ ) is the part of G,.r = 0; Eq. (|85|) 
is a combination of the parts of Geg = and G^e = 0, and the part of Grr = 0. We use ([8^) in the form 



fs + fi + fg = 0. The justification for this is that /a, and /g must be zero in a stationary solution |31|. The 
integration constant must be zero therefore, when ( |8^ ) is integrated. Similar arguments justify integration of (|84|). 
With (^9|) -(^5|), and considerable manipulation, we end up with 

Qo = ^ (364'(i)5t*(i) - 19Af5t*(i)9t*(i)) (86) 
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(87) 



S2 ^ -Ldt-^^^^dtt-^^^^ 



126 



(88) 



Notice that Qi and Q2 do not vanish, and hence diverges at r ^ cxd. It is important to understand that this 
is not incompatible with with asymptotic flatness. The value of G^"^^ must fall off as in AF coordinates, so Q*'^^ 



(2) 

can diverge if there is a compensating divergence in ^'r^- There must, in fact, be a divergence of this order, since 
the source term in ( |4^ ) turns out to diverge. In practice, numerical computations with divergent quantities are to 
be avoided. It is useful, therefore, to exploit the fact that the second order wave function is not unique. (See the 

discussion at the end of Sec.ffl) We now introduce an alternative second order wave function 4*^3 j by 



RW 



(2) 
rad 



(2) 
RW 



(89) 



where 



1 



'—rad 



2016 V 



The wave equation for ^'j. 



(2) 
ad' 



■^(*rad) — '^RW + ^(^rad) — S^d , 



(90) 



(91) 



has a source term that is well behaved at r ^ 00, and in AF coordinates, we then have from (^) and (|86|)- (|9^), that 



12rG(^)=vI/(2) +0(r-i). 



and thus 



(92) 



dPower _ 1 / d¥-^^ e^^^rad' 
dn ~ 10247r y dt ^ Y dt 

Integration over all angles then gives us the total power 

1 / e^dM?. 



Power 



ad 



480 \ dt 2 dt 



(93) 



(94) 



VI. SUMMARY AND DISCUSSION 



For convenience, we repeat and summarize here the main results of the paper. The wavefunction 



6 (2r + 3M) 



2(r - 2M)(i/f ^ - rdrK^^y) - 2(r - 3M)K^'^ + 6 IrK^''^ + ^L.^Ml^r^O^.G^^) _ 2/i«) 



(95) 



is Moncrief's p3| wavefunction. It is constructed completely from first order perturbations of the 3-geometry on a 
hypersurface, and its time derivative can be found from the first order perturbations of the 3-geometry and extrinsic 
curvature of the hypersurface. This wave function is invariant with respect to first order gauge transformations, and 
satisfies the Zerilli equation 



92 



Q^2 Qj.2 



V[t) 



= 



(96) 
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where V{r) is given in (E3). From (EJ), (|89|) , (|90[) , our second order wave function is 



(2) 
rad 



6 (2r + 3M) 



(97) 



|^i^(2) + (!_?M)(r2 5^^,(2) _ 2^(2)) I _ (l44vl/(i)ai\l.(i) - 76MatVl;(i)a,v|/(i) 



(98) 



where the exphcit form of Q^yf is given in ( [4l|) . The wave function ^J:^^ gives the second order equivalent of the two 



(2) 



important advantages of the Moncrief function 'f^^^: (i) it is gauge invariant, as spelled out in Sec. [II and (ii) its 
value and its time derivative can be found directly from the first and second order perturbations of the 3-geometry 
and extrinsic curvature of a hypersurface, as described in Sec. IV. The second order wave function satisfies a Zerilli 
equation 



dri 



y{r) = 'Srad . (99) 



The source term Sra.d is well behaved at infinity, and is given by (H), (|0l) and The outgoing solution for ^f^^ 

has the asymptotic behavior of a function of retarded time t — r* , and the gravitational wave power carried in the 
perturbations is given by the simple prescription in (p3|). 

The scheme outlined in the above paragraph gives a formalism that is definitive and complete (for the case of even 
parity, axisymmetric, quadrupole perturbations). This formalism involves only gauge invariant variables, and requires 
the solution of a wave equation with a source that is well behaved at spatial infinity. Due to the gauge invariant 
nature of the wave function, the Cauchy data for the wave equation can be constructed immediately from an initial 
value solution given in any gauge. 

The method given here constitues a reformulation of the fixed-gauge approach used by Gleiser et al. [p^-p^ . The 
gauge invariant approach appears to offer advantages both in the organization of computations, and in the conceptual 
clarity of the gauge invariant variables. 
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